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Abstract 

In a static environment, optional participation and cooperator agglomeration 
are found to be beneficial for the occurrence and maintenance of cooperation. 
In the optional public goods game, the rock-scissors-paper cycles of different 
strategies yield the oscillatory persistence of cooperation but not different 
levels of stable cooperation. In this paper, by incorporating population den- 
sity and individual mobility into the spatial optional public goods game, we 
study the coevolutionary dynamics of strategy updating and benefit-seeking 
movement. With low population density and slow movement, an optimal 
level of cooperation is easy to be reached. Both the increase of population 
density and the speed-up of free-floating of competing agents will lead to the 
decrease of cooperation. A log-log relation between the levels of cooperation 
and the free-floating probability is found. Theoretical analysis indicates that 
the decrease of cooperator frequency in the present model should result from 
the increase of the interactions between different agents, which may originate 
from the enlargement of the average size of the interaction groups or the rise 
of the speed of random- movement. 
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1. Introduction 

Social dilemmas describe a conflict situation existing between a rational 
individual maximizing its own benefit and a social group pursuing collec- 
tive wellbeing[U |2]. For example, as a sheepherder enjoys herding in the 
common grass the joint effort to attain a sustainable environment suffers a 
risk. The gas emission of a factory ensuring its prosperity makes the green- 
house problem become more serious, which in turn do harm to the further 
development of the factory. As everyone faces the temptation to exploit the 
public goods and make no contribution to it, the immense benefit that can 
only be got by all the people's cooperation becomes unattainable. This poses 
a challenging problem about how the individual selfishness can lead to the 
occurrence and maintenance of cooperation commonly found in reality? To 
answer this question, various mechanisms have been introduced by scientists 
since Darwin [31 Hj. 

Most studies on the evolution of cooperation among selfish individuals 
are based on game models, the Prisoner's Dilemma game(PDG)[5l El 13, [8], 
the snowdrift game (SG)P, QUI [HI U21 H3], and the public goods game 
(PGG)[H1 |15l [16], to name just a few. The iterated PDG models trans- 
actions between two agents, in which one's contribution favors the other but 
not itself. Although the total income would be the highest if they both co- 
operate, an agent tends to defect to maximize its own profit. Therefore, the 
Nash equilibrium is to defect at all times. The PGG is an extension of the 
PD to an arbitrary number of agents. In the original PGG, a group of indi- 
viduals have the opportunity to make an investment into a common pool or 
not. An even distribution of returns irrespective of one's contribution makes 
the defectors have an advantage over the cooperators and defection becomes 
a dominating strategy in such an evolutionary system. To refrain from get- 
ting stuck in the deadlocks of mutual defection, a third strategy, termed a 
loner's strategy [T71 [TBI EH ED], has been introduced into the original public 
goods game. 

In the public goods game with loners, also termed optional public goods 
game (OPGG), the agents have the option to participate in the public goods 
game or not. Those who join the public goods game get a cooperator's payoff 
Pc or a defector's payoff Pjj, and those who do not join the public goods 
game get a loner's payoff P^. Because a loner only gets a small but fixed 
payoff, it can win over a group of defectors but will be defeated by a group of 
cooperators. Therefore, an endless rock-scissors-paper cycle occurs: Loners 
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will invade a population of defectors with fewer cooperators. Cooperators 
will thrive in a population of loners. A population of cooperators will be 
intruded into by defectors. 

In the optional public goods game, although the existence of loners help 
the cooperators keep from being doomed, the oscillation of the levels of the 
three strategies indicates that a higher and stable level of cooperation is 
difficult to be reached in such a system. But in the real world, cooper- 
ation often becomes a dominant strategy in animal and human activities 
and an environmental change will lead to the occurrence of different lev- 
els of cooperation. The environmental conditions include the structured 
space |22l 1231 122 ESI EE], the population density [23 123 EH1 EOJ , and the 
mobility of the individuals [HH [32], [331 H2]. To find out the mechanisms de- 
termining the occurrence of different levels of cooperation in the real world, 
it is a real need to generalize the model presented by Hauert et al. and 
incorporate the environmental conditions into the original OPGG [T7] . 

In natural and human society, the structured space is not full-occupied 
and random and purposive movements often occur [351 EEl 1371 ? ]. Such as 
the migration of birds, the floating of a boat and the movement of a train. 
Z.H.Liu has studied the influence of population density and individual mo- 
bility on epidemic spreading [39]. C.P.Roca et al and J.Y.Wakano et al have 
studied the roles of mobility in the improvement of cooperation in the PGG 
and in the ecological PGG respectively [101 ED]- Motivated by the work done 
in partially occupied rigid networks, in this paper, we incorporate population 
density and individual mobility into the original OPGG introduced in ref.|17j 
and play the game in a square lattice with the Moore neighborhood. The 
main findings of the study are as follows: 

(1) With a predefined free-floating probability £ , the cooperator fre- 
quency fc is determined by population density p. There exists a transition 
point ptr, below which fc increases with the rise of p while above which fc 
decreases with the rise of p. With a predefined p, fc is determined by £o- 
Increasing £ leads to a monotonic decrease of fc- 

(2) Considering the size distribution P(S) of individual components, we 
find that the power-law relation between fc and £o is related to the occurrence 
of the largest component. Before the largest component S max ~ iV occurs, as 
p increases, P(S) changes from an exponent to a power-law and the slop of fc 
as a function of £o decreases with the rise of p. After the largest component 
Smax ~ N occurs, as p increases, P(S) changes little and the slop of fc as a 
function of £o has no change with the rise of p. 
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(3) Taking an eye of the decrease of cooperation in the present model, 
the effect of increasing £o with a fixed p is similar to the effect of increasing 
p with a fixed £o- The theoretical analysis shows that the change of the 
frequencies of different strategies in the present model should be determined 
by the possible collisions of the agents with different strategies. Both the 
increase of p and the increase of £ would result in more such collisions. The 
more collisions of the agents, the lower cooperative levels of the system. 

The paper is proceeded as follows. In Section 2, we introduce the op- 
tional public goods game with purposive and random movements in a spatial 
setting. In Section 3, simulation results about the evolution of competing 
strategies and the state of local agglomeration are presented and the rela- 
tionship between them is discussed. In Section 4, the extinction thresholds 
are analyzed with mean field theory, and the possible relations between the 
levels of cooperation and the individual collisions are described theoretically. 
Section 5 summarizes the paper and gives an outlook of future studies. 

2. The model 

A population of N agents is distributed over a square lattice with side 
length L and the Moore neighborhood (i.e., degree k — 8), each agent on 
each site. For N < L 2 , the population density is defined as p = -p. For 
p < 1, an agent can move around and occupy the firstly found empty site. In 
the present model, there exist two co-evolutionary processes: the change of 
personal strategies and the movement of the agents. Once the initial position 
and the adopted strategy of each agent are set, the system can evolve as 
follows. 

In the evolution of personal strategies. Initially, there exist three kinds 
of agents: cooperators (C), defectors (D) and loners (L). At each time step, 
firstly, each agent interacts with its neighbors and get a payoff Pq (for co- 
operators), Pd (for defectors) or Pi (for loners). Owing to the neighboring 
restriction, the interaction group in the present model should be n < 9. As- 
suming n = n c + Ud + n L , in which nc, no, n^ represent the number of 
cooperators, defectors and loners in the interaction group respectively, the 
payoffs for the agents with different strategies are 



rn c 



-1, 



(1) 



Pd = 



n c + n D 
rn c 



(2) 



n c + n D ' 



4 



Pl = a, (3) 

in which r (> 1) is the multiplication rate and a generally satisfies < a < 
r — 1, which keeps a loner's payoff to be lower than the payoff of the agents in 
a group of cooperators and higher than the payoff of the agents in a group of 
defectors. If nc + no = 1, the cooperator or defector will get a loner's payoff. 
After all the agents have attained their payoffs, they will make decisions on 
whether they should update their strategies or not. In the updating process, 
an agent i compares its payoff with a randomly chosen neighbor j's and adopt 
j's strategy with probability 

= i + QiPi-Pj+r)/^ ( 4 ) 

in which r and k represent the cost of strategy change and the environmental 
noise respectively. 

The above updating mechanism shows that, in an interaction group, the 
evolution of the mixed strategy profile is determined by the scores of differ- 
ent strategies. Even if the group size is constant, the number of competitive 
agents, i.e. cooperators and defectors, will vary with time. Such that the 
oscillations of the population composition will occur in a fully-occupied net- 
work setting. 

In relation to the mobility of individuals. For a loner, because it only 
relies on a small but fixed payoff and gets nothing from a benefit-seeking 
competition, it is not so attractive to stay at its current site and at each time 
step it will move randomly with probability For cooperators and defec- 
tors, they make their decisions on whether they should leave their current 
sits or not based on their satisfaction to the present situation and the possi- 
ble free-floating, which are described by the purposive-moving probability £ g 
and the free-floating probability £o respectively. The value of £o is predefined 
and the value of £ 9 is determined by the wealth of an agent, which is an 
accumulated payoff in the latest T time steps. For an agent i, its wealth is 

t 

t'=t-T+l 

For a cooperator or a defector, if its wealth w is less than 0, with probability 

&7D = & + fo, (6) 
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Figure 1: Frequencies of cooperators (circles), defectors (squares), and loners (triangles) 
as a function of multiplication rate r in a square lattice with size M = L X L = 200 x 200 
and Moore neighborhood. Other parameters are a = 1, r = k = 0.1, £l = 0.1, £o = 0.01, 
El = 0.001 and (a)p = 0.3, (b)p = 0.5, (c)p = 0.8. The results are obtained by averaging 
over 10 runs and 1000 times after 10000 relaxation times in each run. 

in which £ fl = 1 — e w , it leaves its current site and finds an empty site to stay 
on. 

Therefore, in the present model, there are two mechanisms which de- 
termine the evolution of cooperation. The payoff determines the replicator 
dynamics while the wealth determines the mobility dynamics. The great 
concern of us is how the evolutionary dynamics in a three-strategy game are 
affected by the population density and the free-floating probability. In the 
present model, to refrain from the random disappearance of loners, at ev- 
ery time step, the competitive agents will occasionally become loners with a 
small random-flapping probability e^. 

3. Results and discussions 

In this section, we will focus on the roles of population density and free- 
floating probability of cooperators and defectors in the change of cooperation. 
Following the work done in ref. [T7], in Monte Carlo simulations we choose 
the loner's payoff a — 1, the cost of strategy change r = 0.1 and the envi- 
ronmental noise k = 0.1. Throughout the paper, the following parameters 
are also predefined: the size of the square lattice M = L x L = 200 x 200, 
the random-moving probability of loners £l = 0.1 and the random-flapping 
C -> L (or D ->■ L) probability e L = 0.001. 

Figure 1 shows the frequencies of cooperators, defectors and loners as a 
function of r in a square lattice with different population densities p = 0.3, 
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Figure 2: Frequencies of cooperators (circles), defectors (squares), and loners (triangles) 
as a function of multiplication rate r in a square lattice with p = 0.3 and (a)£o = 0.01, 
(b)£o = 0.05, (c)£o = 0.1. All the other parameters are the same as those in fig.l. 

0.5 and 0.8. As what has been found in a fully-occupied regular network[25j, 
there exist two extinction thresholds rp and tl- For r < r^, loners perform 
better than cooperators and defectors so that all the agents become loners in 
the final steady state. For r D < r < r L , the three strategies coexist. The rise 
of p is found to be beneficial for defectors and loners but not for cooperators. 
For r > Tl, the survival of cooperators is favored and the loners go into 
extinction. 

In all the three cases with different p, the values of the extinction threshold 
7*d are the same, = cr + 1 = 2, which is also the same as that in a 
fully-occupied regular network|41J. Such a result comes from the fact that, 
with a small multiplication rate r < 2, even if all the competing agents are 
cooperators, the payoff of them, Pc = r — 1, is less than a loner's payoff 
a — 1. Therefore, the rock-scissors-paper cycle does not occur and the 
system is stuck in the state where all the agents are loners. However, as we 
consider the extinction threshold Tl, we find it increases with the rise of p. 
For p = 0.3, tl ~ 4. For p = 0.5, Tl ~ 5.5. For p = 0.8, ~ 6.5. 

The rise of population density can effectively leads to the rise of the 
probabilities of immediate interactions between two agents. Such an increase 
is more possible to occur in the system where the agents can move freely. 
In fig. 2, we give the frequencies of different strateg function of r 

with predefined population density p = 0.3 and different £o- Comparing the 
results in fig. 2 with the results in fig.l, we find that, as to the change of 
the frequencies of different strategies, the effect of increasing £ f° r fixed p 
is similar to the effect of increasing p for fixed £0(7^ 0). As £ increases, the 
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Figure 3: Frequencies of (a) cooperators, (b) defectors, and (c) loners as a function of pop- 
ulation density p for r — 2.5 and different £ = 0.01 (circles), 0.05(squares), 0.1 (triangles). 
All the other parameters are the same as those in fig.l. 



cooperator frequency fc decreases while the defector frequency fu and the 
loner frequency fi increase within a large range of r. As £o increases from 
0.01 through 0.05 to 0.1, the extinction threshold ri increases from 4 through 
5 to 5.5 accordingly. 

The results in fig.l and fig. 2 indicate that, the increase of the meeting 
probability, whether it comes from the rise of population density or the rise 
of free-floating probability, will lead to the decrease of cooperation. 

It is instructive to ask: can we draw a conclusion that the fewer the 
chances to meet, the higher the levels of cooperation in the present model? 
To examine whether an optimal level of population density exists, in fig. 3 
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Figure 4: Time-dependent frequencies of cooperators with r=2.5, £l = 0.1, El = 0.001 
and (a) £ = 0.01, p = 0.07(black), O.lO(red), 0.25(blue); (b)f = 0.01, p = 0.30(black), 
0.35(red), 0.60(blue); (c) p = 0.3, Co = O.Ol(black), 0.05(red), O.lO(blue). 

we give the frequencies of the three strateg ICS cLS cL function of p for different 
£o- With a small and fixed £o = 0.01, there exists a transition point pt r . For 
p < p tr , the defectors do not exist in the population. The cooperators increase 
while the loners decrease with the rise of p. For p > p tr , the cooperators 
decrease while the loners increase with the rise of p. Within this range, the 
defectors firstly increase and then decrease with the rise of p. Increasing 
the free-floating probability leads to the decrease of the levels of cooperation 
but not the disappearance of the transition point. As £o ranges from 0.01 
to 0.1, the transition point p tr changes from 0.28 to 0.1 and the cooperator 
frequency at the the transition point decreases from 1 to 0.86 accordingly. 
The existence of the transition point indicates that, there exists an optimal 
level of population density, with which the system will reach the highest level 
of cooperation. 

To have a close eye of the evolution of cooperation below and above 
the transition point p tr , in fig. 4(a) and (b) we plot the time-dependent 
frequencies of cooperators for different p. For comparison, we also plot the 
time-dependent cooperator frequencies for different £o in fig. 4(c). For p < 
Ptr, as the time passes fc becomes stable. The change of p only leads to the 
change of the average value of cooperator frequencies but not the fluctuations 
of it. For p > p tr , the change of p not only results in the change of fc but 
also the stableness of cooperation. Increasing p leads to large fluctuations 
of the levels of cooperation. Figure 4(c) shows increasing £ also leads to 
large fluctuations of the levels of cooperation. Such results indicate that, 
whether the decrease of cooperator frequency results from the increase of p 
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Figure 5: The size distribution of individual components in the evolved system with r=2.5, 
£l = 0.1, £l = 0.001 and different population density p — 0.3(circles), 0.4(squares), 
0.5(triangles). 
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Figure 6: The size of the largest component as a function of population density p in the 
evolved system with r=2.5, £l = 0.1 and sl = 0.001. 



or the increase of £ , the large fluctuations of the strategies are detrimental 
to cooperation. 

To find the relationship between the frequencies of strategies and the 
population patterns, in fig. 5 we plot the size distribution of individual com- 
ponents for different p. As p increases from 0.3 through 0.4 to 0.5, the size 
distribution changes from an exponential distribution through a power-law 
distribution to a single quite large component accompanied by a power-law 
distribution. 

In fig. 6 we give the size of the largest component as a function of p. As p 
ranges from 0.1 to 0.3, the size of the largest component has little change with 
the rise of p. As p ranges from 0.3 to 0.5, the size of the largest component 
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Figure 7: Log-log plot of the cooperator frequency fc vs £ with P = 0.3, 0.4, 0.5, 0.6, 0.7. 
All the other parameters are r=2.5, £l = 0.1 and sl — 0.001 . The fitted lines satisfy the 
equation fc = a^Q b . 

has a sharp increase with the rise of p. For p > 0.5, nearly all the agents are 
in the same component and the largest component has little change with the 
rise of p. 

Comparing the results in fig.3 with the results in fig. 5 and fig.6 we find 
that, for low population density, the agents are scattered and the levels of 
cooperation are somewhat high. For high population density, nearly all the 
agents are in the same component and the the levels of cooperation become 
lower. Therefore, the results in fig.3 can be understood as follows: For 
quite low population density, the leave of the unsatisfactory agents makes it 
become impossible for the defectors to exploit the cooperators and the group 
of defectors will finally be doomed by the random-moving loners. And for the 
cooperators who leave the competing group, because of the low population 
density, it is not easy for them to find another cooperator to collaborate 
with and they will finally become solitary agents whose payoffs are equal 
to the payoff of a loner. The occasional C — > L (or D — > L) random- 
flapping can make the loners refrain from extinction. Therefore, under such 
circumstances, only cooperators and loners are left in the final steady state. 
But for high population density or a high free-floating probability, it is easy 
for defectors to intrude into cooperator clusters and the average level of 
cooperation will accordingly decrease with the rise of p or the rise of £o- 

To have a deep understanding of the roles of free-floating in the increase or 
decrease of cooperation, in fig. 7 we plot fc vs £o for different p. From fig. 7 we 
find that, for low population density, the frequency of cooperators decreases 
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sharply with the rise of £ - The rise of p makes such a changing tendency 
become ease. As we fit curves to data points in fig. 7, it is found that fc 
and £o satisfy the equation fc = a^Q b , in which a=0. 28769, 0.37352, 0.39333, 
0.37861, 0.34281 and 6=0.23855, 0.12845, 0.087778, 0.07006, 0.069364 re- 
spectively. 

Comparing the results in fig.7 with the results in fig. 6, we find that the 
slope of the fitted line in fig.7 is closely related to the occurrence of the largest 
component. As the size of the largest component S max is quite small, the 
slope of the fitted line is steep, which indicates that the change of £o can 
quite effectively affect f c as the agents are scattered. As S rnax has a sharp 
increase, the slope of the fitted line quickly becomes gentle. As S rnax ~ N, 
the change of p has little effect on the change of S max , and accordingly the 
slope of the fitted line no longer changes with the rise of p. Such results 
indicate that, in the present model, the value of b in equation fc = a^Q b 
contains the information of evolutionary patterns. We may effectively figure 
out the size distribution of the components from the slope of the fitted line 

of/o(&). 

The above simulation results indicate that, as to the decrease of cooper- 
ation, the role of increasing £o is similar to the effect of increasing p. The 
change of cooperation in the present model should come from the change of 
the meeting probability of the agents. Both increasing £ and increasing p 
can effectively increase the meeting probability of the agents with different 
strategies, which makes the cooperators easy to be exploited by the defectors 
and accordingly leads to the decrease of cooperation. 

4. Theoretical analysis 

4-1- Relationship between extinction threshold tl and average size n of in- 
teraction groups 

In the present model, due to the purposive movement and free-floating of 
the agents, the interaction partners in the competing group vary with time, 
which will lead to the similar results in well-mixed populations where all the 
agents are possible to be chosen as interaction partners in the competing 
process. In the following, we make a mean field analysis of the replicator 
dynamics and give an approximation of the extinction threshold in the 
present model. 

According to the payoff function, in a randomly chosen group, a defector 
always gets a higher payoff than a cooperator. But in mean field analysis, the 
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system does not evolve according to such payoffs but the averaged payoffs 
for cooperators or defectors which are obtained by averaging over all groups. 
In the present model, not all the interaction groups have the same size. In 
theoretical analysis, we take the average size of the interaction groups as the 
interaction group size, which satisfies n = 9p. 

Suppose in the well-mixed population, the frequencies of cooperators, de- 
fectors and loners are fc, Jd and ft respectively, which satisfy the condition 
fc + Id + /i = 1- As that in ref.[T7j, the average payoffs for defectors, 
cooperators and loners are 



_ f 1 fn 

Pc = P D -{r-l)fr 1 + r -\^±-l, (8) 

Pl = o. (9) 

As the multiplication rate r increases to the extinction threshold r — Tl, the 
loners become extinct and we get /l = and fc + /d — 1- In such 
the payoffs for defectors and cooperators become 

tT — 1 

Pd = ~^r L f C} (10) 
n 

n — 1 „ 1 . . 

Pc= ^^r L f c +-r L -l. 11 
n n 

On the condition that Pc = Pd, where cooperators and defectors can coexist, 
we obtain the extinction threshold 



r L = n. (12) 

According to the above equation, the extinction threshold in the 
present model is only related to the average size of interaction groups. The 
rise of population density will lead to the rise of the average group size and 
thereafter the increase of the extinction threshold. For example, as popu- 
lation density p increases from p = 0.3 to p = 0.8, we can estimate that 
tl should increase from 2.7 to 7.2. Compared with the simulation results in 
fig.l, it is found that the theoretical analysis is in accordance with simulation 
data. 
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4-2. Relationship between free-floating probability £ and levels of cooperation 

In simulation results, it is observed that the free-floating of the agents has 
great impact on the change of cooperation. In the following, by theoretical 
analysis, we will give a picture of what may be the possible reasons for the 
occurrence of such an impact. 

In the OPGG, the frequencies of cooperators, defectors and loners are 
determined by the payoffs of the strategies. In well-mixed cases, the payoff 
for each agent is determined by the group size and the status of the agents 
in the same group. The larger the group size, the more possible for two 
agents to meet each other at a time. Just like that in well-mixed cases, in 
a mobile environment, although all the agents are on the sites of a network, 
they are possible to meet each other within a period of time. The faster 
the free-floating of the agents, the more possibly the agents meet each other 
within the period of time. Therefore, from the view point of the meeting 
probability of two agents, the effect of increasing the speed of free-floating is 
the same as the effect of increasing the group size. In the following, we will 
firstly give a functional relation between the group size and the free-floating 
probability, and then give a theoretical analysis of the relationship between 
the free-floating probability £ and levels of cooperation. 

In the present model, because all the agents are arranged on the sites 
of a rigid network with Moore neighborhood, the number of agents in the 
same group should be less than or equal to 9. Considering the effects of 
increasing the group size and the speed of free-floating on the increase of 
meeting probability, we define the boundary conditions of £o and n as: for 
£o ~ 0, n ~ 9p; for £o ~ 1, n ~ 9. Therefore, the following functional relation 
between n and £o is adopted, 



Just like that in a well-mixed case[T7l HTJ, for predefined p and £ > the 
difference of the payoff between a defector and a cooperator satisfies the 
equation 




(13) 




(14) 



Pd = Pl = cr. Therefore, the above equation becomes 
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,n(p,£ )-l r 1 - /. 



L 



1 + -^)^7T^°- (15) 

For < fi < 1, the solution of the above equation is the same as the solution 
of the following equation 



[l-h) + (r-l)(l-f L )f n L (Mo) ' 



n(p, Co) 



0. 



(16) 



Suppose 



WlXMo)) = (l-/L)+(r-l)(l-/ i )/ L " ( ^ 0) " 1 -- 7 -yT(l-/L ( ^° ) )- (17) 
For the three cases of f L = 0,f L = l and / L — > 1, we obtain 

F(0,n(p,e )) = l-^ TT , (18) 

F(l,n(p,£ ))=0, (19) 



F(f L ,n(p,^)) « (1 - r)(n(p,&) - 1)(1 - /l) 2 . (20) 

From the equations of (18), (19) and (20) we find that, within the range of 
r > 2 and n(p, £ ) > r > as /j, changes from to 1, the value of n(p, £ )) 
changes from F(f L ,n(p,£ )) > to F(f L ,n(p,£ )) < 0. As we check the 
sign of 9 F Qf2 ,n ^ , w e find that it changes at most once within the range of 
< /l < 1. Such results indicate that, within the range of < /l < 1, there 
exists a single solution /£ r for equation (16). For fi > /| r , the payoff for 
cooperators is less than the payoff for defectors. For f L < f l L r , the payoff for 
cooperators is greater than the payoff for defectors. 

The above theoretical analysis indicates that, in the present model, in- 
creasing £ will lead to the increase of fi and accordingly the decrease of fc, 
which is in accordance with the simulation results found in Fig. 3. 
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5. Summary 

In the face of the choice to win everything or get nothing in a public 
goods game, it is a reasonable choice for the individuals to drop out from the 
game and enjoy a small but fixed gain, which yields the rock-scissors-paper 
cycles of different strategies in the optional public goods game. The oscilla- 
tory persistence of cooperation in the optional public goods game displays 
the sustain-ability of cooperation, but it does not tell us on which conditions 
the system will evolve to slow oscillations and different levels of stable coop- 
eration can be reached. By incorporating population density and individual 
mobility into the original OPGG, it is found that both the stableness and 
the improvement of cooperation are connected to the degree of crowdedness 
and the speed of free-floating of competitive agents. 

For low population density and slow free-floating of competitive agents, 
the departure of unsatisfactory agents from competing groups make it easy 
for scattered cooperators to agglomerate into independent cooperator clusters 
which results in the expansion of cooperation. The agents who are left on the 
original sites are more possible to form a single defector cluster, which will 
finally be doomed by the random-moving loners. For low population density 
and fast free-floating of competitive agents, the exchange of neighbors often 
takes place and the defectors are easy to invade the independent clusters of 
cooperators. Numerous small components where cooperators and defectors 
coexist occur and the levels of cooperation decrease. For high population 
density, because of the percolation effect, nearly all the agents are merged 
into the same component. The defectors uniformly expand into cooperator 
territory and the evolutionary dynamics in the original OPGG is recovered. 
The relationship between the levels of cooperation and the free-floating prob- 
ability is found. 

The present model indicates that the existence of different levels of coop- 
eration in the real world should be related to the environmental conditions, 
including population density and individual mobility. In the future, simi- 
lar environmental conditions should be considered in the game models with 
continuous strategy variables and the generalization of these conditions is a 
special interest of ours. 
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